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Abstract: Theories of massive gravity inevitably include an auxiliary reference metric. 
Generically, they also contain an inconsistency known as the Boulware-Deser ghost. Re- 
cently, a family of non-linear massive gravity actions, formulated with a flat reference 
metric, were proposed and shown to be ghost free at the complete non-linear level. In 
this paper we consider these non-linear massive gravity actions but now formulated with a 
general reference metric. We extend the proof of the absence of the Boulware-Deser ghost 
to this case. The analysis is carried out in the ADM formalism at the complete non-linear 
level. We show that in these models there always exists a Hamiltonian constraint which, 
with an associated secondary constraint, eliminates the ghost. This result considerably 
extends the range of known consistent non-linear massive gravity theories. In addition, 
these theories can also be used to describe a massive spin-2 field in an arbitrary, fixed 
gravitational background. We also discuss the positivity of the Hamiltonian. 
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A. 4 The Boulware-Deser analysis 



1. Introduction and summary 

Generically, theories of massive gravity are plagued by the Boulware-Deser (BD) ghost 
instability at the non-linear level |IJ gj. Recently, significant progress has been made 
towards constructing massive gravity theories that avoid this instability. In addition to the 
metric g^ u , theories of massive gravity inevitably include another rank-2 symmetric tensor 
fpu, henceforth called the reference metric. This is due to the fact that the interaction terms 
that can be formed from the metric alone, tr g = 4 and det g, cannot be used to construct 
a mass term. Most of the recent work has focused on the case of a flat reference metric, 
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essentially, fnv = flaw I n particular, a two-parameter family of actions was proposed in 
||, H for this case by demanding the absence of the BD ghost in what's known as the 
decoupling limit. One of these actions was demonstrated to be ghost-free more generally 
at fourth order in perturbation theory in jjj. The full two-parameter family of actions were 
then shown to be free of the BD ghost instability at the complete non-linear level in || 
based on the reformulation given in |H|. For complementary work see ||, [|. 

In this paper we consider non-linear massive gravity actions constructed with a general 
f^p and extend the proof of the absence of the BD ghost given in || to this case. This 
generalization is motivated by several considerations. First, there is no reason to insist that 
a theory of massive gravity always refer to a flat reference metric. For example, one may also 
consider dS or AdS metrics. Second, forcing f^ u to be flat constrains the classical solutions 
of the metric g^. For example, non- Minkowski homogeneous and isotropic spacetimes 
were argued to be excluded in |i~Q| . A possible resolution of this problem is to allow 
for more general f^ u |J. A third motivation is that, from a theoretical standpoint, it is 
more satisfying to promote f^ u to a dynamical field with its own kinetic term than to 
have a "frozen-in" reference metric. The resulting theory would resemble the bi-metric 
construction of [11, 12]. For a dynamical fn V to be consistent, it is important to first verify 
that the mass term which was ghost free for flat remains so for a general f^ u . 

It should be emphasized that although the discussion in this paper is formulated in 
the context of massive gravity, the analysis applies equally well to generic massive spin-2 
fields. For example, g^ u could also represent a neutral massive spin-2 meson in a fixed 
gravitational background f^. 

At the linear level, massive gravity theories with a few simple non-flat /^'s have 
already been considered. The linear Fierz-Pauli theory of massive gravity jD^, [l4j in a flat 
background has been generalized to linear massive gravities in de Sitter and anti de Sitter 



spacetimes |15|, [16|, [L7|, [T^] , and to FRW backgrounds |19|, |2(], |21|] . In these constructions 
f^v plays the role of the background FRW metric. Our work explicitly shows that it is 
possible to construct non-linear extensions of these theories that are free from the BD ghost 
instability even for a general f^ v . Namely, the non-linear ghost-free massive gravity actions 
proposed in [|| remain ghost-free when constructed with respect to a general f^ v . 

In this work we consider the massive actions of § as reformulated and extended 
to general f^ v in ||. In this reformulation, the two-parameter family is regarded as an 
extension of a simpler, "minimal" massive action. Each free parameter is associated with a 
higher level of non-linear complexity. The simplicity of the minimal model is instrumental 
in constructing the proof of the absence of the BD instability. Moreover, for this model 
the constraint equations can be solved explicitly, making it possible to study issues such 
as the positivity of the Hamiltonian. Once the proof of the absence of the BD ghost is 
constructed for the minimal model, we find that the exact same construction holds for the 
more complicated two-parameter family of actions. 

Our analysis is based on the ADM formulation of gravity [22] . In the ADM language, 
the BD ghost is a consequence of the absence of the Hamiltonian constraint. We show 
that in the models considered here, such a constraint exists. With an associated secondary 
constraint (see 22]), this is enough to eliminate the BD ghost mode. 
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The paper is organized as follows: Section 2 starts with a review of non-linear massive 
gravity with a general reference metric fn U - We discuss the Boulware-Deser ghost problem 
in non-linear massive gravity and present precise criteria for avoiding it. We then review 
the specific two-parameter family of actions considered in this paper. In section 3 we show 
the absence of the BD ghost in the minimal massive action by obtaining the Hamiltonian 
constraint and arguing for the existence of an associated secondary constraint. We then 
discuss the positivity of the Hamiltonian. In section 4, the proof of the absence of the 
BD ghost is extended to the complete two-parameter family of massive actions and the 
Hamiltonian constraint is determined. The results are briefly discussed in section 5. In the 
appendix we review the ghost issue in linear and non-linear massive gravities, including 
the original analysis of Boulware and Deser |]J . 

2. Review of massive gravity and the Boulware-Deser ghost problem 

In this section we discuss the ghost problem in massive gravity, reviewing the Boulware- 
Deser argument and the caveat by which it can be avoided. We also review the two- 
parameter family of potentially ghost-free massive gravity actions formulated with respect 
to a general f^ u . 

2.1 General structure of non-linear massive gravity 

A generic covariant massive gravity action for the metric is obtained by adding a 
non-derivative potential term V(g^ 1 f) to the Einstein-Hilbert action Q, 



S m = Ml J d 4 x^g- [R(g) - m 2 V{g- l f)] . (2.1) 



fnv is a non-dynamical rank 2 tensor that is needed to construct generally covariant, non- 
derivative functions of the metric. The coupling of the metric g^ to matter is taken to be 
the same as in GR in order to preserve the weak equivalence principle. Below, we will have 
more to say about the role of f^. 



Such a generic non-linear massive gravity action (2.1) typically contains a ghost, i.e., a 
physical mode with negative kinetic energy which in the quantum theory results in negative 
probability states. The origin of the ghost is easy to understand (see the appendix for a 
detailed review of the ghost problem). In general relativity, a scalar component of the 
metric is potentially a ghost, but is eliminated by the equations of motion. The addition 
of a potential energy term to the GR action generally results in this component becoming 
an independent dynamical degree of freedom. Then the theory has six propagating modes: 
the five polarizations of the massive spin-2 graviton and the ghost. 

At the linear level, the ghost problem is avoided by the mass term proposed by Fierz 



and Pauli |Tj|, 14]. A necessary (but not sufficient) requirement for the action fl2,l|) to be 
ghost-free is that, when expanded to quadratic order in metric fluctuations huv = g^u — g^, 
the potential V should reproduce the Fierz-Pauli (FP) mass term in the background g^, 
provided one takes g^ y — f/ivi 

Sfp = - X -Mlm 2 f d'x^g [h$ h» - {h»f ] , (2.2) 



- 3 - 



where, hy = g^ p h pv . At the linear level, this particular choice of relative coefficients 
between the terms decouples the ghost by making it infinitely massive. However, it was 
shown by Boulware and Deser [l], §1 that the ghost sixth mode generally reappears at the 
non-linear level, leading to speculation that a non-linear theory of massive gravity may not 
exist. We will discuss this further below. 

For the purpose of identifying the massive excitations, f^ u must be equated to a back- 
ground metric, f^ u = g^ v . Then the background field equations reduce to the GR equations 
with a shifted cosmological constant and f^ is a solution for a given source, say, T^ u . For 
this reason f^ v is often referred to as a "background metric". But, given /, the non-linear 
theory will also have classical solutions in which g differs appreciably from / in some re- 



gions of spacetime, see for example, |25, [27|. Any such solution can be regarded as 
a background g with fluctuations h' pu around it, although the action for these fluctuations 
may no longer have the Fierz-Pauli form ( |2.2[) . Thus at the non-linear level, one could 
consider fluctuations around background metrics other than f^. For this reason we refer 
to f^ u as the "reference metric" , rather than a background metric. The physical metric of 
spacetime is still g^ u . 

Of the ten components of f^, four are gauge degrees of freedom, removable by gauge 
fixing general coordinate transformations. This is made explicit in the parameterization, 

hv ~dx^ ab dx»- ( ' 

The 4> a are interpreted as Stiickelberg fields or as Goldstone modes associated with the 
breaking of general covariance |28|| . The remaining six components contained in / are 
non-dynamical. Possible choices for f pu are: 

Flat reference metric: Most of the recent work on massive gravity has focused on 
ffiu = Vfiv- In the unitary gauge this gives f pv = rj^. For this choice, (|2.2D is the original 
ghost free Fierz-Pauli mass term [|1|, |l4|] for metric fluctuations around flat spacetime. 
Later, the generic instability of the non-linear theory (|2.1[) was shown by Boulware and 
Deser Q, § for this case, although their analysis also applies to general /. The actions 
recently proposed in [||, || also belong to this class, where the fields <ft a ( |2.3| ) played an 
important role in the construction. The absence of the Boulware-Deser instability at the 
complete non-linear level was proved for these actions in ||. 

FRW reference metric: The quadratic action (|2.2[) is also known to be free of the 



Boulware-Deser ghost instability when f„ v is a de Sitter or anti de Sitter metric [15, 16 



17, 18 1, or more generally, an FRW [19, 20, ^] metric 1 . However, consistent non-linear 



extensions of such quadratic actions had so far remained undetermined. 

General non- dynamical reference metric: In this paper we demonstrate the consistency 
of the non-linear massive actions proposed in §, when extended to general f^ v ||. Such 
an extension is not only natural, but is also necessary to obtain a larger and potentially 
more viable class of solutions. 



It turns out that, in some regions of parameter space, these theories may suffer from instabilities quite 
distinct from the Boulware-Deser problem, even at the linear level. However, these do not necessarily reflect 
an inconsistency of the theory jl^, [HJ . 
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Dynamical reference metric: It is appealing to complete the theory by including dy- 
namics for fp V [11, 12 1 . That this can be done consistently in the context of bi-metric 



theories of gravity will be demonstrated in an accompanying work 
2.2 The ADM formulation of general relativity 

The physical content of gravity and its propagating modes are easily identified in the ADM 
formulation p2j which is based on a 3 + 1 decomposition of the metric, 



N = (-9 00 r 1/2 , Ni = g oi , //,, . (2.4) 

The N and iVj are the lapse and shift functions respectively. In this parameterization, 

/ -1 iV j \ 

9^ = N- 2 [ . . , (2.5) 

where, N J ' = 7 J iV^ and 7^7,^ = d~ l k . Denoting the momentum canonically conjugate to 
jij by 7T 1 -?, the Einstein-Hilbert action in terms of these variables becomes (we ignore all 
boundary terms in what follows), 

S = Ml j d 4 x [n^dtjij + NR° + N,R l ] . (2.6) 
The are functions of 7^ and it 13 but are independent of the iV„ = (N,Ni), 



R° = ydeT 



R i = 2 v / det 



7 V^ 



( 7T*- 



\/det 7 



(2.7) 

The six components of 7^ are potentially propagating modes in the sense that their equa- 
tions of motion obtained from ( |2.6D , as well as those for their conjugate momenta tt^ , 
involve time derivatives (so that the Euler-Lagrange equations for 7^ are second order 
in time). Since a single propagating mode involves a field component and its canonically 
conjugate momentum, the six potentially propagating modes are described by the 12 func- 
tions (7y , 7r ,,J : ) . However, in the theory defined by ( |2.6[ ) not all of these are independent. 
To see this, note that the appear linearly as Lagrange multipliers, hence their equations 
of motion are four constraints (the "Hamiltonian" and "momentum" constraints) on the 
remaining fields, 

i?o(7,vr) = 0, i^(7,7r) = 0. (2.8) 

These constraints can be used along with the four general coordinate transformations to 
eliminate eight of the 12 functions, in favor of two remaining pairs. These pairs are the two 
propagating modes of GR, describing the two polarization states of the massless graviton 
at the non-linear level. In particular, the scalar ghost is not part of the physical spectrum 



m 



Of the 12 equations of motion for (7^, tv 13 ), four reduce to Bianchi identities while 
another four determine the N^. The remaining four equations describe the propagating 
modes. 
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2.3 The Boulware-Deser ghost 

Boulware and Deser used the ADM formalism to study the physical content of non-linear 
massive gravity and argued that, generically, massive gravity has six propagating modes. 
The sixth mode is the ghost that was avoided in the linear FP action but reappears at 
the non- linear level They also found that the non- linear theory had a pathological 

non-positive Hamiltonian. Let us summarize their analysis here. 

In ADM parameterization, the massive gravity action ( |2.1| ) becomes, 

S m = Ml j d A x [TT ij d m + NR° + Nitf - mV( 7 , N, N h fa)] , (2.9) 

where V = V det 7 NV and coordinate transformations have been used to set = / M „ 
( ]2.3D . Since V is a non-linear function of the , the lapse and shift are no longer Lagrange 
multipliers and their equations of motion, 

dV 

R»( 1 ^)=m 2 V»{ 1 ,N,N i J), with Fe — , (2.10) 

no longer constrain 7^ and n 13 . Instead, these equations can be solved for the Nn in terms 
of (jij , 7r y ) . After eliminating the Nn in this way, one is left with twelve equations for the 
twelve dynamical variables (7^ , ) , hence the theory contains six propagating modes. In 
particular the sixth, ghost mode that was avoided in the linear FP theory has re-emerged 
as a propagating mode. Boulware and Deser Q argued that, since in massive gravity V' 
is always non-linear in the N^, the sixth mode cannot be avoided. 

As an explicit example of a non-linear mass term, considered the FP mass fl2.2j ) 
with = 7]^ u , where now h^ v = g^ u — r]^ u is no longer treated as a small fluctuation. 
This analysis is reviewed in the appendix. They found that in the linear approximation 
the Hamiltonian constraint eliminates the ghost, hence the linear FP theory is indeed 
consistent. However, at the non-linear level, for the mass term considered, there is no such 
constraint. Thus, to reiterate, the Boulware-Deser ghost instability is due to the loss of 
the Hamiltonian constraint at the non-linear level. Finally, in [||, || it was concluded that, 

• The massive theory has six rather than five degrees of freedom, and hence contains 
a ghost. 

• The Hamiltonian of the massive theory is not positive definite. 

• In the limit m — > the Hamiltonian diverges, hence this limit does not exist. 

There is, however, a caveat in the arguments of Boulware and Deser. As pointed out 
in 0], avoiding the Boulware-Deser instability does not, in fact, strictly require linearity of 
the theory in the N^. Rather, it is enough that one combination of the four Nn equations 
of motion fl2.10|) becomes a constraint on the (jij,^). Based on this observation, we 
formulate the criteria for the absence of the ghost in the following subsection. 
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2.4 Criteria for the absence of the Boulware-Deser ghost 

The caveat in the Boulware-Deser argument can be stated as follows. While the potential 
V(g^ 1 f) = V(N, .ZVj,7, /) is a non-linear function of the iV M , suppose there exist potentials 
V for which the equations of motion depend only on three combinations of the N^, 

n r = n r (7V,7V',7) , r = 1,2,3. (2.11) 

That is, the equations ( 2.10|) take the generic form, 

#*(7,tt) = m 2 V^( 7 ,n r ). (2.12) 

Then, in principle, three of these equations can be used to determine the n r in terms of 
(7, 7r). Substituting the result into the remaining equation gives a constraint on the (7, n) 
that may have the right form to eliminate the ghost field. Finally, one also needs a second 
constraint to eliminate the variable canonically conjugate to the ghost field. 

The linear Fierz-Pauli theory is linear in the lapse but not in the shift. Thus the N 
equation of motion provides a modified Hamiltonian constraint while the N equations are 
not constraints, but determine the iVj in terms of the (7, ir). It is natural to expect that this 
feature extends to a ghost-free non-linear theory, especially since the ghost is a scalar. Thus 
the functions rii(N, Nj, 7) can be regarded as the counterparts of N in the massive theory, 
consistent with the 3-dimensional general covariance maintained in the ADM formulation. 
Assuming that the functional relationship is invertible (as it should be), one can determine 
the Ni as functions of the nf N(N,nj, r y). Then the massive action (|2.9| ) can be expressed 
in terms of the combinations rii, 

S[N,Ni) = S[N, nj {N,N i} >y)] . (2.13) 

Now consider the equations of motion, 



SS _ 5S 
SNi Srij 



drij 5S _ 5S 

n6N~~ ' SN = SN 



SS 
Snj 



Sf = 0, (2.14) 



n 5N 



where the subscript on the vertical bar indicates the variable held fixed in the process of 
variation. This leads to the equivalent equations, 



6S_ 



— 

n SN 



0, (2.15) 



which are linear combinations of the iV„ equations. Based on these equations one can now 
formulate a nested set of criteria for the existence of a Hamiltonian constraint. 

1. As described above, for a constraint to exist, the n« equations of motion should 
depend on N^ only through the three combinations nj. Then they can be used to 
determine the rij in terms of (7^ , tt 13 ) . 

2. The equation also must involve only the and be independent of N so that, 
given the n % solution, it becomes a constraint on (7^,71"*''). For this to be the case, 
the massive action in the form S, i.e., when regarded as a functional of N and rii, 
must be linear in N. 
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3. The action S also contains the term N^R 1 where iVj = Ni(N,iij,j). Linearity of S in 
N then implies that the expression for N{ in terms of nj must be linear in N. 

Note that the minimal coupling of the metric to matter is linear in the lapse and shift. 
If this were not the case, the constraints of even massless GR would be violated. Thus 
criteria specified here are not modified by the presence of the minimal matter coupling. 

In this paper we show that, for the massive gravity theories described in the next 
subsection, once requirements 2 and 3 are satisfied, then 1 follows automatically. This 
guarantees the existence of a Hamiltonian constraint associated with the N equation of 
motion. We also argue for the existence of a non-trivial secondary constraint (subsequently 
proven in |23f |) as, simply, the non- linear extension of the known secondary constraint in 
the linear FP theory (see, e.g., |30|). These two constraints eliminate the canonical pair 
corresponding to the Boulware-Deser ghost, reducing the number of propagating modes 
from six to five. 



2.5 Non-linear massive gravity actions with general /, 



(JLV 



In principle, the above criteria might be used to construct a theory of massive gravity that 
does not suffer from the Boulware-Deser instability. In practice, this was not so straight- 
forward. 2 Potentially ghost free actions were first constructed for f^ u = r/^ u following a 



very different perturbative argument. It was observed in |^] that the 4> a in (2J3) are the 
Goldstone bosons associated with the breaking of general covariance by the mass term. 
Then for f^ u = 7]^ u and g^ u = rj^ + h^ u , an analogy with the Goldstone- vector boson 
equivalence theorem in gauge theory implies that, in the high energy limit, the dynamics 
of massive gravity is mirrored in the dynamics of the Goldstone sector, particularly, in the 
"longitudinal" mode of the <j) a fluctuations. Thus the ghost of massive gravity appears as a 
ghost in this longitudinal mode. Being a scalar field in flat spacetime, this is a much easier 



setup to investigate. One may attempt to constrain V using this correspondence p2 |. 

The breakthrough came with the work of de Rham and Gabadadze Q and de Rham, 
Gabadadze and Tolley [0] who used this approach to construct a two-parameter family of 
massive actions for f^ u = rj^. The two free parameters are denoted by and 04. These 
actions were shown to be ghost free in this high energy limit, the "decoupling limit". To 
go beyond the decoupling limit, Q carried out an ADM analysis of the = 0:4 = 
model to quartic order in the metric perturbations and demonstrated the absence of the 
Boulware-Deser ghost to that order. 

The presentation of the actions given in Q is convenient for a perturbative analysis. 
However, the expressions that multiply the parameters a n contain mixtures of terms with 
different levels of non-linear complexity. To demonstrate the absence of the Boulware- 
Deser ghost at the full non-linear level it is helpful to use the reformulation of || in which 
different levels of non-linearity are disentangled. Using this reformulation, the absence of 
the Boulware-Deser ghost at the non- linear level was proven for = r\^ v in ||. 



2 In hindsight, these criteria are powerful enough that they can determine the complete form of the 
non-linear action as will be discussed elsewhere |3l|. 
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In this paper we extend the ghost analysis to any general, non-dynamical f^ v using the 
presentation of massive gravity actions given in H. Let us briefly review this formulation. 
The basic building block of non-linear massive gravity is the square-root matrix y g~ 1 f 
[|, where y g^ 1 f\/ g~ l f = g^fxv The terms appearing in the massive action are 
identified as elementary symmetric polynomials of the eigenvalues of this matrix. They 
sum up to a "deformed determinant" (for related ideas see [^, ^7|). The antisymmetry 
property of this structure allows one to generalize the reference metric from flat to any f^ v 
and still remain within the same 2-parameter family of actions. 

The simplest non- linear massive action with zero cosmological constant is given by @ , 



M 2 

±v±p 



j d 4 x^g R-2m 2 (tr yV 1 / - 3 ) • ( 2 - 16 ) 



We will refer to this as the minimal massive action. The most general non-linear massive 
action can be written as 

3 

(2.17) 

n=0 

The efe(X) are elementary symmetric polynomials of the eigenvalues of X. For a generic 
4x4 matrix they are given by, 



r r 3 , 

5 = M 2 / d A x^~g R + 2m 2 Y, Pn e„( V^F?) 

J L „ — n 



e (X) = 1, 

ei(X) = [X], 

e 2 (X) = i([X] 2 -[X 2 ]), (2.18) 

es(X) = i([X] 3 -3[X][X 2 ]+2[X 3 ]), 

e 4 (X) = -^([X] 4 - 6[X] 2 [X 2 ] + 3[X 2 ] 2 + 8[X][X 3 ] - 6[X 4 ]) , 

e fc (X) = for k > 4, 

where the square brackets denote the trace. Of the four f3 n , two combinations are related 
to the mass and the cosmological constant, while the remaining two combinations are free 
parameters. Setting the cosmological constant to zero and the parameter m as the mass, 
the four f3 n are parameterized in terms of the «3 and 04 of Q as (for n = 0, . . . ,4), 

Pn = (-1)™ Q(4 - n)(3 - n) - (4 - n)a 3 + a 4 ) , (2.19) 

The minimal action corresponds to /3 2 = ^3 = supplemented by /3o = 3, j3\ = — 1 to get 
a zero cosmological constant contribution from the potential term. We will start with this 
minimal theory in the analysis that follows. 

3. Absence of the BD ghost in the minimal massive action 

In this section we show that the minimal non-linear massive gravity action ( |2.16| ) satisfies 



the criteria outlined in section and thus there exists a Hamiltonian constraint on the 
dynamical variables. In addition, we argue for the existence of an associated secondary 
constraint. Thus this action does not suffer from the Boulware-Deser ghost instability. We 
solve the constraints explicitly and discuss the positivity of the Hamiltonian. 



-9- 



3.1 Enforcing the criteria for the existence of the Hamiltonian constraint 

In the ADM formulation the Lagrangian for the minimal massive action ( [2.16 ) becomes, 

C mi n = 7r ij d aij + NR° + N'Ri - 2m 2 ^/d^N ftr - 3) , (3.1) 



where, using the parameterization (2.5) for g^ v , one gets, 

N 2-i f= ( ~foo + N l f lo -hj + N'fu \ 

Since the action contains the square root of this matrix, it is highly non- linear in the N^. 
Hence it could potentially propagate a ghost sixth mode according to the Boulware-Deser 
argument. But, as discussed in section 2A, this can be avoided if the four equations of 
motion happen to depend only on three combinations of the lapse and shift, say n l (Nu), 
leaving a single constraint to eliminate the sixth, ghost mode. 

We show now that this is indeed the case for the minimal action fl3.l|) . First, we 
identify the appropriate functions n\ This is achieved by imposing the criteria for the 
absence of ghost discussed in section 2A. In fact, we will only have to impose criteria 2 
and 3. Then 1 follows automatically. 

Criterion 3 requires that the expression for iVj in terms of the rij be linear in N, 

N* = j + Nd l . (3.3) 

The c l and d l are functions of ri 1 and jij but are independent of N. They will be determined 
in what follows by demanding that the action, when written in terms of nj and N, must 
be linear in N. Using (|3.3[), (3.2) takes the form 



AT 2 5 - 1 / = E + iVEi + iV 2 E2 (3.4) 
where the matrices Eo, Ei and E2 are independent of N. To write them compactly, define, 

a M = -fotM + c l fip . (3.5) 



a i \ I ° 1 



/ d l fio d l fu \ 

Ei = , ■ , • • ( 3 - 7 ) 

V ~(d l fioc l + a d l ) -(dd l f tj + d» 0i ) I 




Criterion 2 of section |2.4| requires that the mass term, when written in terms of N 
and Hi, must be linear in N. For the minimal massive action, this is satisfied if the matrix 
\J g~ 1 f has the following form 3 , 



Ny/g-^f = A + NE , (3.8) 



3 Note that this is more restrictive than requiring the linearity of the tr{N^/ g~ x f) in N, but leads to 
simple systematics that satisfy requirement 1 automatically. 
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where the matrices A and B are independent of N. Demanding that this expression fl3.£| ) 
be consistent with N 2 g~ l f as given by (3.4), determines A and B as well as c 1 and d\ 
Explicitly, comparing ( |3.8| ) and (|3.4| ) gives, 



A 



En 



E 9 



and AB + BA = Ei . 



(3.9) 



Let us consider the first two equalities in (p7g|) . Using ( p\qj it is easy to verify that these 
imply, 



1 



a Q 

-anc* 



-c a,- 





D^Cf-^fio 







(3.10) 



Here 3 /j ? - = /j,- and we have introduced, 



x = an — c'a/ 



(3.11) 



The expression for A = VEq follows since Eo is a projection operator, Eq = xEn, hence 
A = E.o/y/x. In the expression for B, the square-root matrix D is defined with an extra 
factor of \fx for later convenience. 

Before proceeding further, we note a very important property of the matrix D. Ac- 
cording to ( 3. lip it has the form D = \J S 3 f where both S and 3 / are symmetric matrices. 
By rewriting D as yl + (S 3 f — 1) and then expanding in powers of (S 3 / — 1), it becomes 
obvious that ( 3 fD) T = 3 fD. In terms of components, this means 



fikD k j — fjkD k i . 



(3.12) 



This identity will be used often in the following analysis. 

Now let us consider the third equality in fl3.9|) . Using ( 3.10| ), one can compute AB + BA 
and compare the result with Ei in (3^7). Using ( [3.5| ) and the property of D given in ( 3,12| ), 
one obtains the following relation between c* and d l , 



d 1 = D l 



c k - CrTfio 



(3.13) 



Guided by the case of flat f^v considered in M, we introduce the variables n l so that 4 , 



n k = c k 



(TYfio- 



Then, (|3.13j) reduces to, 



d 1 = D\ n k . 



Substituting for d l in ( |3.11 ) gives a matrix equation for D, 

y^D = ^( 7 -i - Dnn T D T ) 3 f , 



(3.14) 
(3.15) 

(3.16) 



4 This choice simplifies some of the equations, but is not unique. For instance, we could have also chosen 
n 1 = c\ For a different choice, see section 3.5 below. 
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which will be solved below in terms of n\ Thus d, d l and D l ■ can be determined entirely 
in terms of the n l and jij . This proves that indeed there exist modified shift variables n l 
in terms of which NyJ g~ 1 f is linear in N and given by ( |3,8[ ). 



In the proof that follows, we need only the condition (3.16), and not its explicit solution. 
However, it is important that this solution exists and can be used to show that the relation 
(3.23) is invertible. Thus we take a moment to derive the solution. Squaring both sides of 
(3.16) and moving the D-dependent terms to one side gives, after using ( p. 12 ), 

!)',(/,!>', "/,/,•• with, Q l j =x5 l j + n l n m f mj . (3.17) 

On multiplying both sides by Q, this becomes (D Q) 2 = (7 _1 A f)Q- Taking the square root 
and rearranging gives, 

D = (y/r^Q)Q~ 1 - (3.18) 
The inverse matrix Q _1 is easily obtained by noting that (nn T 3 /) 2 = {n T 3 fn)nn T 3 f. Then 
one finds 

Q- 1 = -(1- M~ 2 nn Tz f) , (3.19) 

x 

where M 2 = —foo + fok( 3 f~ l ) kl fio is the lapse function of f^ u . Equations (|3.18| ) and ( 3.19j ) 
give the explicit solution for D in terms of n\ 

Before moving on, let us point out that our final expressions naturally involve the 
ADM parameters of f^ u . In analogy with the ADM parameterization of g^, we define, 

M = (-/ 00 )-V2 } Mi = f 0i , % = Sij. (3.20) 

We also define M l = (^f ) %3 Mj. Then, the variables defined in ( |3.5| ) become, 

a = M 2 + n l Mi , a, = %n l , (3.21) 

and in terms of the lapse M, the x of ( |3.11 ) is simply, 

x = M 2 - n k f k i n l . (3.22) 



To recapitulate, we have identified three variables n l such that the Lagrangian ( |3.1| ) 
written in terms of these variables is linear in iV and hence satisfies criteria 2 and 3 for 



the existence of a Hamiltonian constraint, as outlined in section 2.4. The functions n % are 
related to the lapse and shift variables of g^ u and through 5 , 

AT* = n * + M l + N D\ n k , (3.23) 

where the matrix D is given by ( |3.18j ) and ( 3.19j ) above. Note that n % parameterize the 
difference between the shift functions of g^ u and fav- In the following section we will show 
that, with this choice of n l , criteria 1 of section |2.4| is automatically satisfied. 



For / M „ = rjftv, in a perturbative relation between Ni and rn was used to demonstrate the absence of 
the BD ghost to quartic order in the fluctuations. To compare this perturbative relation with our result we 
expand ( |3,23| ) to quartic order around flat spacetime using ( 3. IS ) . After lowering the indices and rescaling 
rii by 2, one gets, 

Ni = m + i^iVn, - ^SNh/rij + ±(-h/ + | h/h^n, . 

In Q] the third term comes with a coefficient 1/8 and the fourth term is absent. Thus it appears that at 
low orders this relation is not unique. 
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3.2 The Hamiltonian constraint in the minimal massive theory 



Now we consider the N and n % equations of motion ( 2,15 ) and show that they do not depend 
on N. Hence they give rise to a Hamiltonian constraint. Incorporating (|3,8| ), ( 3,10| ), ( 3.15| ) 
and ( |3.23|) into the minimal massive theory (3.1), leads to an action in terms of the n % 
which is linear in N, meeting requirements 2 and 3 of section 2.4, 



TT ij d t 7ij + NR° + Ri [n* + M i + N D\ 



n 



-2mVdet7 [Va? + N^tx D - 3N] . (3.24) 



Thus the N equation of motion (2.15) is independent of N by construction. We now show 
that the n % equations are also independent of N as required by criterion 1. To get the n l 
equations, one needs, 

d ,— 1 dn l 



dn k 



-4= "7: 
fx 



^(v^trD) = --ln^ 



(3.25) 
(3.26) 



The first line easily follows from (3.22). In the second line, we have first used ( 3.1C| ) 



to re-express the left hand side in terms of the square root matrix and then S tr yE 



|tr(v-E 5E) to evaluate the derivative. On using (3.12), the right hand side follows. 
Then, varying C m i n in ( [3.24 ) with respect to n k gives, 



Ri + 2m 2 v / det7 



nj fj 



d 
dn k 



in 1 + ND\n l 



(3.27) 



Note that the expression in the square brackets is the Jacobian matrix of ( [3.23D . From 
( 2.14 ) it is then obvious that the expression in the parenthesis is indeed dC m i n / dN l . The 
Jacobian matrix is generically invertible as can be checked, for example, perturbatively by 
using the expression in footnote (|5|). Hence one gets the n % (or N l ) equations of motion, 



\fx Ri + 2m 2 \J det 7 /y n j = . 



(3.: 



which are independent of ./V as advertised. Using ( 3.22| ), these determine ri 1 in terms of 7^ 
and the conjugate momenta tt 1 -? (contained in Ri), 



-M 



n 



^/Am^t 1 + R k {f- 1 ) kl R l 



As a consistency check, note that this implies, 

4m 4 det 7 M 



4m 2 det 7 + R T3 f- 1 R 



> 0. 



(3.29) 



(3.30) 



so that \fx is real on the constraint surface. 
The N equation of motion is, 



R° + RiD)n j - 2mVdet7 y/xD k k -3 



(3.31) 
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The n % that appear explicitly and through x ( 3.22 ) and D ( |3.18 ) can be eliminated using 



the solution (3.29). Thus the ./V equation becomes a constraint on the 12 components of 
jij and tv 13 and reduces the number of canonical variables to 11. This is the Hamiltonian 
constraint of the minimal massive action. 

To see that this constraint has the correct form to eliminate the ghost, we adapt a 

5 T (2), ij 



parameterization of [22] for 7^ in terms of the six functions 7^ T (2), tJ(2), 7 L (1)1 and 



7 J (1), 

7y = 7? T + %nj + djji + • (3.32) 

Here, 7^ = tj(<% — V~ 2 didj)j T and 7$ = jf + dij L . As the notation implies, 7^ T is 
traceless, transverse and "yj is transverse; hence the above counting. -f T is the trace of the 
transverse part of 7^'. The flat space limit indicates that jJj T , jj and j L carry the massive 
spin-2 graviton while "f T describes the ghost. From the analysis of [22] one can see that 7^ 



appears in R° in the right way to be eliminated by ( |3.31 ), in analogy with GR. One more 



constraint is needed to remove the canonically conjugate variable. 

3.3 The secondary constraint and the absence of the BD ghost 

Now we argue that the Hamiltonian constraint gives rise to a secondary constraint (for an 



explicit proof see |23j) so that the 12 dimensional phase space of 7^ and n 13 has only 10 
degrees of freedom, corresponding to the five polarizations of the massive graviton. 

Before proceeding note that after integrating out the shifts N\ the Lagrangian (3.24) 
remains linear in the lapse N, 

Cmin = K ij d tlij - Uofrij, 7r ij , f) +NC( 7ij y\ /) . (3.33) 

From this one can read off a Hamiltonian, ignoring the usual ADM contribution that can 
be expressed as a boundary term, H = J d 3 x(T-Lo — NC). 

A secondary constraint is obtained by demanding that the Hamiltonian constraint 



(3.31), now summarized as C = 0, is independent of time on the constraint surface. In the 



Hamiltonian formulation this condition is expressed in terms of a Poisson bracket, 

jC = {C, H} « (3.34) 

with H as given above. Now, if {C(x), C(y)} 0, this condition becomes an equation for 
N and does not impose any constraint on the 7^ and Tr tJ . If this were true, as argued 
to be the case in | ]33| [, then a second constraint would not exist. However, in |23j it 
has been shown through explicit calculation that {C(x), C(y)} ~ 0. Then the condition 
dC/dt = is independent of N and thus becomes a second constraint on 7^ and it 13 (with 
H = J d 3 xU ), 

C {2) = {C, #o}~0, (3.35) 

provided the expression for C {2) does not vanish identically. That this is the case can 
be easily shown perturbatively. By construction, the Lagrangian ( |3.33| ) reproduces the 
Fierz-Pauli Lagrangian at lowest order in the fields and for /„„ = r]^ u . Hence, 

H ^< P + O(7 3 ,vr 3 ,5/), C^C FP + 0( 7 2 ,7T 2 ,5f). (3.36) 
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Therefore, one should find that, 



C P)- C £f + Otfy>V), (3-37) 

where it is known that is neither identically zero nor equal to C FP (see, for example, 
pOp ). Moreover, as can be seen from the Fierz-Pauli structure, enforcing dC^/dt = 
will result in an equation for N and no further constraints are generated. For details see 
p3[| . To summarize, this proves the existence of a secondary constraint that is needed to 
completely eliminate the propagating Boulware-Deser ghost mode. 

3.4 The positivity of the Hamiltonian and open issues 

In general relativity, the Hamiltonian corresponding to the ADM Lagrangian (|2.6D is ex- 



pressible as a boundary term [22|. The boundary expression for the Hamiltonian can also 
be derived in a more general setup |34], ^] by considering the Gibbons-Hawking boundary 
terms that have been suppressed in ( |2 .6[ ) . This boundary expression for the Hamiltonian is 
independent of the mass term and remains unchanged in the massive theory. However, the 
mass term gives an extra bulk contribution H m to the Hamiltonian due to the reduction 
in the number of constraints. 

In the massive gravity actions that they analyzed, Boulware and Deser found that the 
contribution of the mass term to the Hamiltonian % m was generically not positive and 
moreover, it diverged in the limit m — > (as reviewed in the appendix). The pathological 
behavior of the Hamiltonian was given as a reason for disregarding massive gravity [ID, |j . 

Here we consider the contribution of the mass term to the Hamiltonian in the minimal 
massive action with a general reference metric. This can be easily computed from C m i n 
in ( |3.24| ) upon imposing the Hamiltonian constraint ( 3.31] ) . Using (3.29) and (3.30) this 
contribution becomes, 



n min , m = My 4m, 4 det 7 + R k ^f- l ) kl Ri — R%M\ (3.38) 

where, M and M l are the lapse and shift functions of Note that for f„ u = rj^ u 

(M = 1, M % = 0), which was the case considered by Boulware and Deser, this is clearly 
positive and well behaved in the limit m — > 0, avoiding the pathologies observed in fij. 
The same applies to any reference metric f^ v for which M l = and M > 0. 

In general, when M* ^ 0, the last term in T-L m in,m appears problematic, but at least in 
some cases this is a gauge artifact. For example, even for a flat f^ u = d^(f) a d v 4> b ri a b, if we 
do not choose the physical gauge <p^ = x^, the Hamiltonian is not manifestly positive. But 
clearly, in this case, the problem is a gauge artifact. It seems that in all situations where 
one can choose a gauge with M % = 0, the Hamiltonian remains positive. 

However, while the minimal massive action with a general reference metric is free of 
Boulware-Deser instability, it seems that it may not always have a positive Hamiltonian. 
In fact, given the known instabilities of massive gravity in, say, de Sitter and anti de Sitter 
backgrounds, we do not expect a massive gravity Hamiltonian to be positive for all possible 
f^v in all regions of parameter space. However, precisely because of the RiM 1 structure, the 
theory is potentially linear in the lapse and shift functions of f^, provided the constraints 
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do not introduce a non-linear dependence through the elimination of the ghost mode. This 
opens up the possibility of consistently promoting /„„ to a dynamical variable. This will 
be discussed in detail in [29]. 

3.5 An alternative set of variables 

It is not apparent that the equations of section 3.1 depend in a simple way on the lapse 
M of fn V . However, the solution ( |3.29| ) shows that the n l are linear in M, so that x 
is proportional to M 2 and the matrix D goes as 1/M. This motivates working with an 
alternative set of variables, 

n* = Mn* , D) = D)/M. (3.39) 
which makes the simple dependence of the theory on M manifest even beyond the minimal 



model. In term of these, (3.23) takes the form, 

N l = Mfj + M i + N D\ h k , (3.40) 

and 

x = M 2 x , x = l- h T3 fh (3.41) 
The defining equation for D, and hence the matrix D itself, is independent of M , 



yflb = \fi~f- 1 - Dhh T D T ) a f , (3.42) 

Then it follows that, expressed in terms of hatted variables, the matrix B in ( |3.10| ) is 
independent of M. The matrix A depends on M in a more complicated way through 
ao = M 2 + Mh l Mi, di = Mh l Mi and c l = Mh l + M i . However, the most general massive 
action (|2.17 ) contains A only in the combinations tr(A), tr(AB) and tr(AB 2 ) (see the next 
section). It is easy to verify that all these are linear in M. Thus, on using the hatted 
variables, the action ( 2.17D becomes linear in both iV (see below) and M. 



4. Absence of the BD ghost in the complete 2-parameter massive action 

We now demonstrate the absence of the Boulware-Deser ghost in the full two-parameter 
generalization of the minimal massive theory with a general reference metric We show 
that the Hamiltonian constraint is maintained even in this case. The analysis is more 
involved but it turns out that the variables identified in the minimal massive model can 
be used in these more general cases without any modifications. The argument for the 
existence of the secondary constraint parallels the discussion for the minimal model. The 
details are given in ]23fl. 



4.1 The 2-parameter action in terms of the new variables 



We now show that, using the same n l of the previous section ( 3.23 ), the general massive 



action (2.17) turns out to be linear in ./V thereby satisfying criteria 2 and 3 of section 2.4 



In the ADM parameterization, the general Lagrangian is given by 

C = ^d tlij + NR° + N'Ri + 2m 2 ^N /3„e„(^ T 7 )] . (4.1) 



\n=0 
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To check if this Lagrangian is linear in the lapse N, let us express the potential in terms of 
the matrices A and IB of ( |3.8| ), Since e n ~ (\/~g~\f ) n > h might seem that negative powers 
of ./V should appear in the Lagrangian. However, due to the property tr(A fc ) = (trA) fc , 



these terms cancel amongst themselves. The potential terms (2.1 g|) then give, 



Ne 1 {y r g~ T f) = trA + iVtrl 



Ne 2 {^g- l f) = trAtrB-trAB + ±iV [(tr B) 2 — tr B 2 ] , (4.2) 
Ne^g- l f) = tr AB 2 - tr ABtr B + \ tr A [(trB) 2 - trB 2 ] 
+ ±iV [(tr B) 3 - 3tr BtrB 2 + 2tr B 3 ] . 

These terms are at most linear in the lapse N and thus satisfy criterion 2. The eo term 
just contributes to a cosmological constant while the e\ term was already considered in the 
previous section for the minimal action. So the terms that remain to be investigated are 
e2 and e^. 

4.2 The Hamiltonian constraint in the 2-parameter theory 

We now compute the n* equations of motion and show that they do not depend on the 
lapse N. The RiN 1 term in the action contributes a term Ri J l k to the equations of motion, 
where, 

is the Jacobian matrix of ( [3. 23 ). Since the Jacobian contains the lapse N, the only way for 



the n l equations of motion to be independent of N is if the contribution of the potential 
terms is also proportional to the Jacobian. This must happen separately for the e 2 and 
terms. We show now that this is the case. 

In the following, we employ matrix notation where n is a column vector with elements 
n % and n T is its transpose. To vary with respect to n fc , consider the A terms first. Using 



(|3.10[) , these terms are, 

trA = \fx , 

tr AB = -n T 3 fDn, (4.4) 
trAB 2 = -^n T3 f D 2 n = --^n T3 f (-y- 1 - Dnn T D T ) 3 f; 



n . 



The last expression is written in two equivalent ways using ( 3.16 ). Varying these gives 



dn k ^ ^ \fx ^ f dn k ' 

-A trAB = -n T 3 f D & - n T 3 f 

an K J an K J an h 



_^ tr AB 2 = (n T 3 fD* n) n T3 f-^- 2^ n T 3 f D (4.5) 

T 3f t,2^ \ „T 'if dn r, /- T3 f n 2 dn 



(n^/^n) n-/ - 2^n-^ 2 
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From the last equality we derive the relation, 



1 



y/xD + —Dnn T3 f 



d(Dn) 
dn k 



y/xD 2 + — l(n T3 fD 2 n) 

x 



dn 
dn k 



(4.6) 



This relation is very useful. We need all derivatives in d(Ne%)/dn k to appear in the 
combination ( [4 . 3| ) . However on direct substitution we will find some d(Dn)/dn k terms 
without the factor N. Equation ( |4.6D allows us to re-express these in terms of d(n)/dn k . 
Now consider the B terms. Using ( |3.10| ) these are, 



trB = v^trD, tr B 2 = ^fx 1 tr D 2 

The variations can be written as, 



tr B 3 = yfx tr£> 3 



(4.7) 



AtrB 

an" 



J_ „T 3 f d(Dn) 
dn k 



dn k 

Combining all these results gives 



j^tiE 2 = -2n T3 fD^, 

on K J on K 1 

^trB3 = -3^n^^l. 



(4.8) 



,9 
dn 1 



Netiy/f^f) = -^=n T3 f^(n + NDn), 
v x 



dn 



|^iVe 2 (^/^ T 7) = n T3 f [D-ttvD] £,( n + NDn), 



(4.9) 



dn 



^Nesiy^f) = -V^n T3 f [D 2 - Dtr D + \t [(tr D) 2 - tr(D 2 )]] ^( n + NDn) 



where (4.6) was used to simplify the last expression. Note that the right hand sides are 
proportional to J l k = dN l /dn k ( [4.3| ) which was a requirement for the n l equations of 
motion to be independent of N. So, finally, varying the general action ( [4.1[) with respect 
to n % gives the TV-independent equations of motion, 



.r 



Pi 4 + h (SjD m m - D 1 



+ ft Vx- 2 {\5\{D m m D n n - D m n D n m ) + W m D™ - D\D m m ] j = . (4.10) 

In principle, these equations can be solved to determine the n % in terms of R{ at least 
perturbatively, although unlike the minimal massive model, an explicit non-linear solution 
may be difficult to obtain. 

The N equation of motion is, 

R° + /»',/)>' + 2m 2 ^[/? + p X yftD\ + \foyfx~ 2 (D\D 3 ) - D^D^) 

+ i/? 3 ^ 3 (D\D^D k k - 3D\D\D k 3 + 2D^D\D k ^ = . (4.11) 

Eliminating the n* in favour of R 1 converts this into the Hamiltonian constraint on 7^ and 
7T 1 -' . This and its associated secondary constraint are enough to eliminate the ghost. The 
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existence of a non-trivial secondary constraint follows from an argument similar to one for 



the minimal model in the previous section. For the explicit computations see [23]. 

This demonstrates the existence of a two-parameter family of non-linear theories of 
massive gravity with general fni> that do not suffer from the Boulware-Deser ghost insta- 
bility. 

5. Discussions 

In this work we have shown that the recently proposed non-linear massive gravity theories 
do not suffer from the Boulware-Deser ghost instability for an arbitrary non-dynamical 
reference metric. This is a generalization of the work in ||] which showed the absence of 
the BD ghost for a flat reference metric. To reiterate, the appearance of the BD ghost is 
due to the absence of a Hamiltonian constraint. We have shown that the massive actions 



( 2.17 ) contain such a constraint and an associated secondary constraint and hence are free 
from this instability. 

The theory discussed here need not necessarily be interpreted as a theory of massive 
gravity, which may or maynot be consistent with observations. It also has an alternative 
interpretation as a ghost free theory of a massive spin-2 field (say, a meson) in a fixed 
gravitational background given by f^. 

Much of the recent analysis of the ghost issue in massive gravity has relied on the 
Stiickelberg formulation and a flat f^ u . In the decoupling limit, this formulation provides a 
powerful tool for studying the ghost content of the theory because of the Goldstone-vector 
boson equivalence theorem. The recently proposed massive gravity actions, which we have 
shown to be ghost-free even for a general f^ v , were constructed by demanding the absence 
of the BD ghost in the Stiickelberg formulation and in the decoupling limit alone || |4j. 

However, away from the decoupling limit, the equivalence theorem is no longer valid 
and, even for a flat f^ u , the ghost analysis in the Stiickelberg formulation becomes sig- 
nificantly more involved. A full analysis of the constraints and gauge conditions of the 
theory must be performed to obtain the physical spectrum in order to identify the ghost. 
Some recent work which does not take these issues into account has suggested that, in 
the Stiickelberg formulation (or relatedly, using a helicity decomposition), the BD ghost 
inevitably reappears away from the decoupling limit, at higher orders in perturbation the- 



ory |36|, 37, 38]. These results are in contradiction with the conclusions of |j] and its 
generalization in the present work. However, a thorough analysis of the ghost issue in 
the perturbative Stiickelberg framework was performed in |8|, This showed that when 
the constraints are taken into account, the BD ghost is indeed absent and the apparent 
discrepancy between the perturbative Stiickelberg approach and the non-linear analysis of 
|| disappears. 

By generalizing the ghost analysis to general f^, the results of this paper open up 
the possibility of finding new and interesting classical solutions to massive gravity theories. 
Moreover, it is known at the linear level that massive gravity in FRW-type backgrounds 
may contain instabilities that are distinct from the Boulware-Deser ghost. As argued in 



1 19, 20, 21 1, these problems might be avoided in the full non- linear theory in a dynamical 
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process. The actions studied here provide a non-linear setup in which this issue might be 
investigated. Finally, the results of this paper provide a first step towards promoting f^ u to 
a dynamical variable and thus creating a consistent bimetric theory of gravity (see, p9[| ). 

An interesting possibility is the realization of ghost free massive gravity within string 
theory setups. While this is not obvious at the level of fundamental string, there is evidence 
that linear massive gravity in AdS background arises within the AdS/CFT framework 
1 39, 40]. It is interesting to check if AdS/CFT could also reproduce the correct non-linear 



generalizations described in this paper. 
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A. Appendix: Further review of the ghost problem 
A.l Absence of ghost in general relativity 

In field theory, a ghost refers to a physical mode with negative kinetic energy. In the 
quantum theory this results in states with negative probability. When the action is not 
in diagonal form in the fields and particularly in the presence of constraints and gauge 
symmetries, the physical content of the theory may not be directly discernible. To identify 
the ghost in such cases, one has to first determine the physical degrees of freedom with 
canonical kinetic terms. Alternatively, the ghost appears in the 2-point function as a mass 
pole with negative residue 6 . This comes handy when propagators are known. 

That metric fluctuations in a modified theory of gravity could easily contain a ghost 
can be inferred by investigating linearized general relativity. Decomposing the metric 
fluctuations h^ u = — rj fW in terms of its traceless transverse (h^ v ), transverse vector 
(a^-), longitudinal (4>) and scalar (s) parts, 

V> = h U + d {» a v) + + zV/wS (A.l) 

one can easily check that the quadratic Einstein-Hilbert action depends only on the five 
components of h 1 - and the sixth scalar mode s, 



Seh = \m; / ,i i .v 



h^ u □ h-L - -s □ s 

pv a 



(A.2) 



The other modes drop out due to invariance under 5h^ u = —d^ v y Obviously s has 
a negative kinetic term and is potentially a ghost, but in general relativity it does not 



For fields 4> a with propagators G a b interacting with external sources J a , the transition amplitude 
(0, out\0, in) j ~ e l $ J GabJ will be less than 1 if J a excite |0, in) into particle states of positive probability. 
But this could exceed unity in the presence negative probability states, implying a positive exponent, 
i(2iri)Res(JGJ,k° = |fc|) > 0, §. 
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survive as a propagating mode. In more detail, s contributes to the graviton propagator 
zA^ ) as a ghost. Indeed it contributes the second term in the saturated /c-space propagator, 



rpl£V ^ (Q) rp*pa _ rp_Lf_W rp* _\_ _ rp rp* 



- [W 



6 



(A.3) 



which obviously has a negative residue at the k 2 = pole. However, for k 2 = 0, and only 
for this value, the ghost is cancelled by a similar contribution from the first term, coming 
from the helicity zero component of the spin-2 field h^~. This is discussed below. The 
cancellation is peculiar to the massless theory and does not takes place otherwise. 

A. 2 Avoiding the ghost in Fierz-Pauli massive gravity 

Metric fluctuations around flat spacetime are made massive by adding the Fierz-Pauli mass 



term to the linearized EH action [|13|, 14, 41 1, 

S FP = Seh [h] - \M 2 v m 2 J d 4 x [h$ h» - a {h») 2 ] . (A.4) 

where, tit, = Tf" p h pv . In this form, the mass term is not gauge invariant and depends on the 
a M and <fi of ( |A.1| ). It modifies the massless propagator ( |A.3| ) by shifting the mass poles to 
k 2 = —m 2 and k 2 = —m 2 = \^-^ m2 - Then, on using k fl T tJ ' 1 ' = (see for example, p2[|), 



rpflf \(m) rp*pa 



^ rp-Lpi/ rp*A. _ _ rp rp* 



k 2 + m 2 ^ u 6 k 2 + m 2 



(A.5) 



Then, as shown below, there is no way of cancelling the wrong sign contribution from the 
scalar part against the spin-2 part and all 6 modes (including the ghost) contribute to the 



propagation. The only way out is to tune a = 1 so that m 2 = oo. This decouples the ghost 
keeping only the 5 healthy polarizations of the massive graviton. 

It was pointed out by Boulware and Deser |l]] that this method of avoiding the ghost 
cannot be easily implemented beyond the linear theory. Then, generically the theory will 
contain 6 propagating modes indicating the reappearance of the 6th mode that was removed 
at the linear level by setting a = 1. We discuss this below. 

A.3 Unitarity analysis of the saturated propagator 



The treatment here follows [43|. In our conventions, a negative residue of the saturated 
propagator T pv Ap Upa T* pa implies the presence of a ghost. Given a 4- vector k^ we can 
construct a set of 4 linearly independent vectors, 

k>* = {k ,k), k" = (-k°,k), e^ = (0,e r ) for (r = l,2) (A.6) 



such that e r ■ e s = S rs and k ■ e r = 0. Note that (kk) = kpk* 1 = (k ) 2 + \k\ 2 > and k 2 = k 2 



In this basis, a generic symmetric tensor can be expanded as 



Tp V {k) = ak»k u + bk»k u + \c rs {e^e u s + e^) + \d(k»k v + k"k*) + \e r {We u r + We?) 

+ y r (k p e u r +k u e?) (A.7) 
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The conservation equation k^T^ = implies b = ak 4 /(kk) 2 , d = —2ak 2 /(kk) as well as 
f r = —k 2 /(kk)e r . The propagator (A. 3) in Einstein-Hilbert theory takes the standard 
from on using T^ v = T^ v — ^(rj^u — ^§-)T. Then in the above parameterization of the 
residue at the zero mass pole is positive, hence the theory is ghost free, 



\TT* 



k 2 =0 



Men - C22I 2 + 2|ci 2 | 2 > 0. 



(A.8) 



In the massive theory, the propagator has the generic momentum space form (A.5). The 
T^~T term is due to spin-2 exchange and gives a positive residue at the mass pole, 



k 2 =—m 2 



ak 2 (l 



k 2 k 2 
(kk) 2 



c 

+ 2 



+ 5IC11 - c 22 | 2 + 2|ci2 1 



+ ^ 1 ! 2 + |e 2 | 2 )^-l) 



k 2 -- 



> (A.9) 



m 2 (- 



+ 



The first line is manifestly positive. In the second line, /c 2 (-^- — l)|i.2__ IT ,2 

J 1 ' K (kk) nk ~ m y m 2 +2\k\ 2 

1) > 0, hence the overall positivity. 

The TT* term in ( |A.5| ) is due to the exchange of the scalar mode s of mass m s which 
is a ghost for any finite mass since, 



6" 



k 2 -- 



ak 2 (l 



k 2 k 2 
(kk) 2 



k 2 =—m 2 . 



< 



(A.10) 



Only for k 2 = this cancels against a contribution in the T ± T ± term resulting in the 
healthy GR expression above. The only other possibility to get rid of this ghost is to take 
m s — , 00, for fixed m. This decouples the scalar ghost and retains only the healthy spin-2 
contribution (|A.9|). This is the Fierz Pauli massive gravity for a = 1. 



A. 4 The Boulware-Deser analysis 

As a specific example, considers a FP-type mass ( pT2| ) where h 
treated as a small perturbation. Then, in ADM variables, 



{IV 



g^u - Vnu is not 



K K ~ (K? = (M 2 " (M 2 + 2/i(l - A^ 2 + NiNrf') - 2NM 



(A.ll) 



where hij = "fij — 5ij and h = ha. This is obviously non-linear in both N and N l . The 
equations of motion for these can be solved to give, 



N = —- 



R° 
m 2 h 



N; 



— [(hi 



-1 



r-l\-l 



(A.12) 



There is no constraint on the remaining variables and the theory contains 6 propagating 
modes, including the ghost. Substituting back in the action one can compute the Hamil- 
tonian density % = n^dfjij — C (ignoring the boundary contribution that is the same as 
in GR), 



H = -[( hij ) 2 -h 2 ] + -hm 2 + -^ 



1 1 

2 m 2 



i? i [(/ i7 - 1 -i- 1 r l ] i ,^' (A.i3) 
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JH noted that the corresponding Hamiltonian is not always positive and that it diverges in 
the limit m — > 0. Hence the conclusion that gravity cannot have a finite range. 

It is instructive to see how the ghost disappears in the linear FP limit. In this case, 
expanding around a flat background, 5N = N — 1, 5N{ = Ni and h{j = 7^ — dij are small 
perturbations and to quadratic order the FP mass term is linear in SN, 

K K - {h^f = (hijf - {hf - Ah5N - 2SNi6Ni (A.14) 

The Ni equations R l (h, 5ir) = m 2 5N l determine the <5iVj. But the 5N equation R°(h, 5ir) = 
m 2 ha is independent of lapse and shift and becomes the modified Hamiltonian constraint. 
The requirement that this constraint is maintained under time evolution, results in a non- 
trivial secondary constraint on the (hij , Sir 11 ) . These 2 constraints reduce the number of 
independent (h, 5ir) components from 12 to 10, implying the existence of only 5 propa- 
gating modes and no ghost mode, consistent with the manifestly covariant analysis of the 
propagator. 
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